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Abstract 

We compute the isospin 1 = and 2 amplitudes for the decay of a kaon into two 
pions by estimating the relevant hadronic matrix elements in the chiral quark model. 
The results are parametrized in terms of the quark and gluon condensates and of the 
constituent quark mass M. The latter is a parameter characteristic of the model, 
that we restrict by matching the results in the two 75-schemes (HV and NDR) of 
dimensional regularization. We find that, for values of these parameters within the 
current determinations, the A J = 1/2 selection rule is well reproduced by means of 
the cumulative effects of short- distance NLO Wilson coefficients, penguin diagrams, 
non-factorizable soft-gluon corrections and meson-loop renormalization. 
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1 Introduction 



The AJ = 1/2 selection rule |T| states that, in the weak non-leptonic decays of 
kaons (as well as of hyperons) the amplitude in which the change in isospin is 3/2 
is very suppressed with respect to that in which the change is 1/2. 

For the decay of a neutral kaon into two pions, the CP-conserving amplitude 
with a final / = state (A/ = 1/2) is measured to be 

Re Ao(K° ^ 2tt) = 3.33 x 10"^ GeV , (1.1) 

and it is approximately 22 times larger than that with the pions in the 1 = 2 state 
(A/ = 3/2): 

Re A2(K° ^ 2n) = 1.50 x 10"^ GeV . (1.2) 

Since a naive estimate of the relevant hadronic matrix elements within the standard 
model leads to amplitudes that are comparable in size, this selection rule has been a 
standing puzzle, the solution of which has attracted a great deal of theoretical work 
over the past 40 years (for a review see, for instance, ref. 0). 

Among the most important steps in the progressive understanding of this selec- 
tion rule, we mention (without any pretense to completeness) the work of Wilson 
Gaillard-Lee and Altarelli-Maiani P who stressed the importance of QCD short- 
distance corrections; the contribution of Vainshtein-Zakharov-Shifman |^ who iden- 
tified in the gluon penguin operators a possible source of enhancement of Aq with 
respect of A2, and the subsequent work of various authors who correctly estimated 
the size of them. 

At the same time, a parallel development took place as Cohen-Manhoar |^ first 
pointed out the relevance for the selection rule of meson-loop corrections in a chiral 
quark model computation, Bardeen-Buras-Gerard [|10| estimated this effect in the 



1/Nc approach with the meson loops regularized by a cutoff and, more recently. 



Kambor-Missimer-Wyler |]rT| in chiral perturbation theory. 
What is then the present status of the rule? 

The two main ingredients entering the theoretical estimate of the selection rule 
are the perturbative determination of the Wilson coefficients of the relevant opera- 
tors and the (non-perturbative) evaluation of the hadronic matrix elements. 
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The short distance part of the analysis can be obtained by a straigthforward 
(albeit challenging) application of renormalization group equation methods. In this 
respect, we have only updated existing analyses by including next-to-leading order 
results [0, 11^ and evaluating the Wilson coefficients at various hadronic scales for 



an updated range of Aq^^, as explained in section 1. 

Most of the uncertainty resides in the second ingredient: the evaluation of 
hadronic matrix elements. The vacuum saturation approximation (VSA) is not 
sufficient — non-factorizable and non-perturbative effects are essential in a dynami- 
cal explanation of the rule — and a model for QCD at low energies is called for in 
order to make a progress in this direction. 

In ref. [|l5| (hereafter referred to as I) we have completed a systematic study of the 
chiral quark model (xQM) [|16], |l^, |l^ and derived the complete O(p^) AS* = 1 chiral 
Lagrangian. The model allows us to determine the coefficient of each term of the 
chiral lagrangian as a function of a few input parameters. All the relevant hadronic 
matrix elements can accordingly be computed and the study of the A/ =1/2 rule 
performed. 

The first ingredient that makes the predictions of the xQM substantially different 
from the VSA results is the inclusion of corrections of order 0{asNc) coming from 
the non-factorizable gluon condensates, first considered by Pich-de Rafael [jl9l (see 
ref. 1 20] for an updating of this analysis). These are truly non-perturbative effects 
which represent a crucial step toward the understanding of the selection rule. Their 
most relevant effect is to give the required suppression of the A2 amplitude. 

To the amplitudes thus obtained, we have consistently added the renormaliza- 
tion induced by meson loops, that we have computed in I. This gives the needed 
long-distance enhancement of Ao while little affecting A2. The meson-loop renor- 
malization also provides the scale dependence of the hadronic matrix elements which 
is necessary in order to make the matching with the Wilson coefficients consistent. 

The idea of matching the scale dependence of the meson renormalization to the 
short-distance /^-dependence of the Wilson coefficients goes back to the work of 



Bardeen-Buras-Gerard |jTO[ (and subsequently |2T[), who, however, applied it in the 
framework of a cut-off dependent regularization and therefore with results that are 
different from ours. 
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All the relevant matrix elements are given in I, to which we refer. We report 
in section 3 those that are most important for the present analysis. In section 4 
we discuss our results at the varying of Aqcd, the matching scale /i and of the 
input parameters: gluon and quark condensates. The xQM leaves us with a free 
parameter: the constituent quark mass M. In order to restrict the possible values it 
can take, we require, as an additional constraint, that the 75-scheme dependence of 
the hadronic matrix elements compensates maximally that of the Wilson coefficients. 

As depicted in fig. 9 of section 5, the combined effects of all contributions thus 
included add up to provide, to a very good approximation, the experimental values 
for the J = and 2 amplitudes. We think that our analysis shows that the A J = 1/2 
selection rule is fully accounted for within the standard model and by our present 
understanding of its non-perturbative aspects. 

This said, a word of caution is perhaps advisable. Various uncertainties related to 
the input parameters are necessarily present in our computation and we will discuss 
them as they occur. In addition, we have to consider the approximations inherent to 
our approach. In particular, higher-order terms 0{p^) in the chiral expansion may 
cause a 20-30% correction [^. This is the systematic uncertainty we ascribe to our 
results. 



2 Effective Quark Lagrangian and NLO Wilson Coefficients 

The quark effective lagrangian at a scale /i < rric can be written as 
G 

^As=i = -^VudV*^J2[^iif^) + ^yiif^)]Qiif^) ■ (2-1) 

The Qi are four-quark operators obtained by integrating out in the standard 
model the vector bosons and the heavy quarks t, b and c. A convenient and by now 
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standard basis includes the following ten quark operators: 

Ql = (Sa'U/3)v_A (%'^a)v-A ' 

Q2 = (sm)v_a M)v-A ' 

(53,5 = (srf)v_AEg(gg)vTA ' , . 

Q7,9 = I (sc?)v-A eg (gg)v±A ' 

«y V±A ; 

where a, [3 denote color indices (a, /3 = 1, . . . , Nc) and Cg are quark charges. Color 
indices for the color singlet operators are omitted. The subscripts (V" ± A) refer to 
7^(1 ±75). We recall that (5i,2 stand for the VT-induced current-current operators, 
Q3_6 for the QCD penguin operators and Qr-io for the electroweak penguin (and 
box) ones. 

The functions and are the Wilson coefficients and Vij the Kobaya- 

shi-Maskawa (KM) matrix elements; r = —VtdVig/VudV*g. The numerical values of 



the Wilson coefficients depend on a^- A recent determination at LEP and SLC 
gives 

a,(mz) = 0.119 ±0.006, (2.3) 
which corresponds to 

Ag)^^ = 350 ± 100 MeV . (2.4) 



The range in eq. (|2.4| ) will be used for our numerical estimate of the amplitudes Aq 
and A2. 

Even though not all the operators in eq. ( p.2| ) are independent, this basis is of 
particular interest for the present numerical analysis because it is that employed for 
the calculation of the Wilson coefficients to the NLO order in as jl^ , [13[| . 

In tables 1 and 2 we give explicitly the Wilson coefficients of the ten operators 
at the scale /i = 1 GeV and = 0.8 GeV, respectively, in the naive dimensional reg- 
ularization (NDR) and 't Hooft-Veltman (HV) 75-schemes. Since Rer = O(10~^), 
the CP-conserving component of Aq and A2 are controlled by the coefficients Zi{fi), 
which do not depend on rrit. 
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a(4) 


250 MeV 


350 McV 


450 McV 




0.113 


0.119 


0.125 


HV 




(0.0320) -0.539 


(0.0339) -0.683 


(0.0355) -0.884 


Z2 


(0.988) 1.30 


(0.987) 1.40 


(0.987) 1.56 


Z3 


0.0060 


0.0108 


0.0210 


Z4, 


-0.0142 


-0.0230 


-0.0389 


Z5 


0.0037 


0.0052 


0.0068 




-0.0128 


-0.0204 


-0.0341 


z-j/a 


-0.0041 


-0.0024 


-0.0012 




0.0087 


0.0140 


0.0234 


zg/a 


0.0016 


0.0068 


0.0140 


zio/a 


-0.0084 


-0.0135 


-0.0223 


NDR 


Zl 


(0.0503) -0.440 


(0.0533) -0.535 


(0.0557) -0.644 


Z-2 


(0.982) 1.23 


(0.981) 1.30 


(0.980) 1.37 


zz 


0.0100 


0.0161 


0.0283 




-0.0281 


-0.0432 


-0.0703 


z^ 


0.0067 


0.0083 


0.0088 


z% 


-0.0281 


-0.0438 


-0.0730 


z-j/a 


0.0039 


0.0122 


0.0218 


--h/h 


0.0120 


0.0208 


0.0378 




0.0102 


0.022.3 


0.0382 




-0.007(3 


-0.011(3 


-0.0180 



Table 1: NLO Wilson coefficients at = 1 GeV in the HV and NDR schemes. The 
corresponding values at = mw are given in parenthesis {a — 1/128). In the HV 
scheme one has 2;3_io (mc) = 0. The coefficients Zi{ij), relevant for the study of CP 
conserving amplitudes, do not depend on rrit. 
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250 MeV 


350 MeV 


450 MeV 




0.113 


0.119 


0.125 


HV 


Zl 


(0.0320) -0.657 


(0.0339) -0.910 


(0.0355) -1.36 


Z2 


(0.988) 1.38 


(0.987) 1.58 


(0.987) 1.96 


Z3 


0.0137 


0.0301 


0.0798 


Z4 


-0.0292 


-0.0540 


-0.115 


Z5 


0.0070 


0.0100 


0.0123 


Z6 


-0.0275 


-0.0515 


-0.112 


z-j/a 


-0.0055 


-0.0030 


-0.0065 




0.0198 


0.0379 


0.0827 




0.0070 


0.0203 


0.0415 


zio/a 


-0.0181 


-0.0330 


-0.0644 


NDR 


Zl 


(0.0503) -0.524 


(0.0533) -0.663 


(0.0557) -0.781 




(0.982) 1.29 


(0.981) 1.39 


(0.980) 1.48 


z-i 


0.0180 


0.0360 


0.0870 


Za 


-0.0471 


-0.0852 


-0.182 


^5 


0.0085 


0.0077 


-0.0129 


z& 


-0.0495 


-0.0947 


-0.226 


Zj/a 


0.0073 


0.0204 


0.0366 


z^joL 


0.0280 


0.0589 


0.143 


Zg/a 


0.0206 


0.0441 


0.0779 


zw/a 


-0.0159 


-0.0267 


-0.0438 



Table 2: Same as in table 1 at = 0.8 GeV. 
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3 The Hadronic Matrix Elements 



The xQM allows us to compute hadronic matrix elements by couphng quarks 
to the chiral Goldstone bosons. The model and the details of the derivation of the 
matrix elements are discussed in I, to which we refer the interested reader. 

The results for the matrix elements 



(g,)o,2 = (2 7r,/ = 0,2|g,|i^0) 



(3.1) 



are given by the equations (8.13)-(8.48) of I, including all contributions of order 
0{N'^), 0{N) and 0{asN). We have used dimensional regularization and, in deahng 
with the 75 matrix, we have worked out the analysis in both the HV and the NDR 
schemes. 

For reference, we report here those elements that are most relevant to the present 
discussion. The contributions to Aq and A2 of the electroweak penguin operators 
Q7-10 are suppressed by the smallness of their Wilson coefficients. 

In the HV scheme we find: 



(Qi)o 

(Ql)2 

(Q2)0 

(^2)2 

(^3)0 
(^4)0 
(Q5>0 



3 



3 



{GG) 



V2 



2 - — (1 - <5(GG>) 



X 



1 + ^ (1 - 5^GG)) 



+ ao(Qi) 

+ a2(Qi) 
+ ao(<32) 

+ a2((32) 



2 m 



{GG) ) + 00(^3) 



Wo = 2^X' + ao{Qe) 



where 



X = y/Sf^ (ml - ml) and X' = X ll - 6 



M2 

X 



(3.2) 

(3.3) 
(3.4) 

(3.5) 

(3.6) 
(3.7) 
(3.8) 

(3.9) 
(3.10) 
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In the NDR: 

(<5i)o 

(Ql>2 
(Q2)0 

{Q2)2 



-X 
3 

V2 



X 



1 + ^ (l - S^OG)) 



-X 



3 



+ ao(Qi) 

+ a2((5i) 
+ ao((52) 

+ a2((32) 



(Q3)0 = ^(X'-5(GG>^)+ao(g3) 

{Qa)o = X' + ao(g4) 

2Mx" + ao(g6) 



(Q6)o ^ 

where 

= X I 1-9 



A2 

X 



(3.11) 

(3.12) 
(3.13) 

(3.14) 

(3.15) 

(3.16) 
(3.17) 

(3.18) 
(3.19) 



The function ao,2(Qi) represent the one- loop mesonic renormahzation of the 
hadronic elements and they are given in I. They are made of polynomial terms, 
generally of the order of 

2 

m 

and logarithmic terms of the order of 



2 2 

m , mi 
In 



m 



2 ' 



(3.21) 



where the masses can be any among m^r, uik and and / is the pion decay 
constant in the tree-level chiral lagrangian. 

The one-loop renormahzation of / is taken into account by replacing / with fi 
in the tree-level amplitudes, which amounts to replacing 1//^ with 1//| multiplied 

by 



1.18. 



(3.22) 
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The corrections of 0{asN) are important. They are parametrized by the value 
of the gluonic condensate: 

AT, {a,GG/n) 

Their most important effect is to reduce the contribution of the operator Qi and Q2 
to the amphtude A2] in fact, 



A2 ^ {zi + Z2) 



(3.24) 



and by taking 5(^gg) — 3, which is what we obtain for the central value of the gluon 
condensate (see eq. ( |3.26| ) below) , the correction is large enough to revert the sign 
of the 1/Nc term and thus suppress the amplitude. 

A relevant contribution to the amplitude Aq arises from the gluonic penguins 
Q5 6, whose matrix elements are directly proportional to the value of (gg) and are 
controlled by M through the suppression factors in ( p.lO| ) and (|3.19|) , that make the 
matrix elements larger for smaller values of M. 

The meson-loop renormalization is sizable. It enhances Aq while little affecting 
y42, as it is needed in order to reproduce the experimental values. An important 
feature of this correction is the additional scale dependence that is introduced in 
the hadronic matrix elements and that matches to a good approximation that of the 
Wilson coefficients. 

For the purpose of comparison with the existing literature, it is useful to intro- 
duce the effective factors 



vVSA ' 



that give the ratio between our hadronic matrix elements and those of the VSA. We 
shall discuss their numerical values in section 4.3. 

3.1 Input Parameters 

The quark and the gluon condensates are two input parameters of our compu- 
tation. As discussed in I, their phenomenological determination is a complicated 
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question (they parametrize the genuine non-perturbative part of the computation) 
and the hterature offers different estimates. 

We identify the condensates entering our computation with those obtained by 
fitting the experimental data by means of the QCD sum rules (QCD-SR) or lattice 
computations. In our discussion we will vary these input parameters within the 
given bounds and obtain a range of values for the amplitudes we are interested in. 

A review of recent determinations of these parameters, together with a justifica- 
tion of the estimated errors, is given in I. Here we only report the ranges that we 
will explore in our numerical analysis. 

For the gluon condensate, we take the scale independent range 

(376 ± 47 MeV)^ , (3.26) 

which encompasses the results of recent QCD-SR analysis . 
For the quark condensate, we consider the range 

- (200 MeVf < {qq) < -(280 MeV)^ (3.27) 



i—GG) 



in order to include the central values and the errors of the QCD-SR |25| and lattice 



estimates 26 



In discussing the scale dependence of our results, it is necessary to include the 
perturbative running of the quark condensate. This can be done in the QCD-SR 
approach by using the renormalizat ion-group running masses m„ -|- m^, the value of 
which is estimated at yU = 1 GeV to be ^7 

mu + md = 12 ± 2.5 MeV (3.28) 

for Ag^D = 300 ± 150 MeV. The error in (|3.28| ) reflects changes in the spectral 
functions. In our numerical estimates, we will take as input values the running 
masses at 1 GeV given by ( |3.28|) . Even though our preferred range of Aq^j-, in eq. 
( p.4|) corresponds to Aq^j-, = 400 ± 100 MeV, and therefore is not that employed in 
ref. 1^^, we feel that we are not making too large an error since the determination 
in eq. (p.28|) is not very sensitive to the choice of Aqcd- 

By taking the value (|3.28|) , we find for the scale-dependent (and normal-ordered) 
condensate 

/>^(1 - 



(qq) (/^) 



(3.29) 
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corresponding ((5,r is a few percent correction) to the numerical values of 

{qq) = -(238 ± 19 MeV)^ (3.30) 



at 1 GeV and 

{qq) = -(222± 19MeV)^ 



(3.31) 



at 0.8 GeV. The error in eqs. (|3.30| )- (|3.31|) corresponds to that in ( |3.28|) . The central 
values in eqs. ( |3.30| )-( p.31| ) are in the lower half of the range ( |3.27| ). 



4 Computing Aq and A2 

The physics entering the AI = 1/2 rule involves very different scales and thus 
different effective theories. As a consequence, any attempt to an explanation must 
necessarily depends on several parameters. 

Our solution has two input parameters, the quark and gluon condensates, as well 
as a free, model-dependent parameter: the constituent quark mass M. In addition, 
the short- distance part depends on several high-energy parameters like the KM 
matrix elements and the masses of the heavy quarks. 

We include all ten operators ( ^.2[ ) even though the effect of the electroweak 
operators is only of a few percents. We proceed step-by-step summarizing our result 
by means of tables and figures. 



4.1 The VSA Approach 

First of all, in order to make possible gauging the progress of our computation 
of the hadronic matrix elements, we report in table 3 the results one would obtain 
by using the VSA for the evaluation of the hadronic matrix elements. 

Several comments are in order. At = 0.8 GeV ~ A^ (that is our preferred 
matching scale) and both in the NDR and the HV schemes, Aq is too small by a 
factor of three while A2 is two times too large. Therefore, the ratio 

uj-^ = Re Aq/ Re A2 ~ 3.7 , (4.1) 

is about a factor six smaller than the experimental value u^'^ = 22.2. 
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/X = 0.8 GeV 


/i = 0.9 GeV 


/i = 1 GeV 




NDR 


HV 


NDR 


HV 


NDR 


HV 




1.11 


0.93 


0.89 


0.77 


0.75 


0.67 


A2 


3.00 


2.76 


3.07 


2.87 


3.13 


2.95 




18% 


15% 


12% 




8% 


7% 


6% 








39%- 


33% 






A^Aa 






4%- 


7% 







Table 3: Matching-scale and 75-scheme dependence of Aq (in units of 10~^ GeV) 
and A2 (in units of 10~^ GeV) in the VSA approach with NLO Wilson coefficients. 
The amplitudes are computed for Aq^j-, = 350 MeV and {qq) given in eq. ( p.29| ). 
The two values quoted for the /^-dependence of the amplitudes correspond to the 
NDR and HV scheme results in the range between 0.8 and 1.0 GeV. 



This estimate is done by using the hadronic matrix element for the penguin 
operators given by the VSA. We have used for the quark condensate the PCAC 
expression given in eq. ( p.29|) . As we shall see in the following, a larger value for Aq 
can be obtained by using larger values for the quark condensate. A popular choice 
is 

\ 1/3 



{qq) (IGeV) 



K 



-260 MeV , 



(4.2) 



^ m,(lGeV) ; 

The literature also offers determinations of the selection rule in which the match- 
ing between Wilson coefficients and hadronic matrix elements is done without meson- 
loop renormalization but at a much lower values of matching (at about 300 MeV, 
for instance |jl9|]). Such a procedure is no longer justifiable in view of the NLO de- 



termination of the Wilson coefficients that shows, for such a low-energy matching, 
the breaking down of the perturbative expansion. 
In table 3 we have defined 

A(0.8GeV) -A,(1.0GeV) 



= 2 



Ai(0.8GeV) +Ai(1.0GeV) 



(4.3) 
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as a direct measure of the scale dependence. The scale dependence of Aq is rather 
large, with A^Aq almost 40% in the NDR scheme. On the other hand, the A2^s 
scale dependence remains within 10%. 

The difference between the HV and NDR results is quantified by 

^NDR _ ^HV 



^NDR _|_ ^HV 



(4.4) 



Contrary to the scale dependence, the 75-scheme dependence is not large for the 
amplitude Ao and A2, remaining below 10% and it is not an issue in this paper. It 
is however very large in other observable quantities like, for instance, e'/e where it 
reaches in the 1/Nc approach the 80% level. 

Because the VSA hadronic matrix elements depend neither on the scaleQ nor the 
75-scheme, the dependences in table 3 are a direct measure of those of the Wilson 
coefficients. 

We now turn to the xQM model determination of the selection rule. 



4.2 Matching and scale dependence 

In matching the short-distance Wilson coefficients to the hadronic matrix ele- 
ments, we combine terms that are scale dependent. Ideally, the scale dependence of 
the Wilson coefficients should compensate against that of the hadronic matrix ele- 
ments to provide a scale independent result. In practice, within our approximations 
and for the central value of Aq^j-,, we find that the scale dependence in our matching 
remains below 20%. In particular, as shown in table 4, the difference in the value 
of Aq between the matching at 1 GeV and that at 0.8 is about 20% in both the HV 
and NDR, to be compared with the values given in table 3, where it is about 30%. 
For A2 it is respectively 12 and 1%. 

As expected, the matching dependence is stronger for larger values of Aqcd- We 
take the central value at 350 MeV as our reference value for the discussion. 

Thanks to the satisfactory scale independence, we can take the matching at any 
values between 0.8 and 1 GeV. Yet it is useful to bear in mind that any number 
quoted in what follows suffers of an intrinsic uncertainty of at least 20% because 

^Except for the perturbative running of the quark condensate. 
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Ag^D = 250 MeV 




^l = 0.8 GeV 


/i = 0.9 GeV 


/i = 1 GeV 




NDR 


HV 


NDR 


HV 


NDR 


HV 




2.31 


2.29 


2.17 


2.15 


2.08 


2.04 


A2 


1.66 


1.56 


1.73 


1.65 


1.80 


1.72 




< 1% 


1% 


2% 




6% 


5% 


5% 




10% - 12% 


\A2 


8% - 10% 


Agio = 350 MeV 




// = 0.8 GeV 


/i = 0.9 GeV 


/i = 1 GeV 




NDR 


HV 


NDR 


HV 


NDR 


HV 


Aq 


2.97 


2.94 


2.66 


2.61 


2.45 


2.39 


A2 


1.60 


1.46 


1.68 


1.56 


1.75 


1.64 




1% 


2% 


2% 


^75^2 


9% 


8% 


6% 




19% - 20% 




9% - 12% 


Ag^D = 450 MeV 




/i = 0.8 GeV 


// = 0.9 GeV 


// = 1 GeV 




NDR 


HV 


NDR 


HV 


NDR 


HV 


^0 


4.26 


4.17 


3.52 


3.41 


3.06 


2.95 


A2 


1.58 


1.34 


1.64 


1.46 


1.71 


1.55 


^75^0 


2% 


3% 


4% 


^75^2 


17% 


12% 


9% 


A;,Ao 


33% - 34% 


A;,A2 


8% - 15% 



Table 4: Same as in table 3 in the xQM approach, for different values of Aq^d. We 
take for the gluon condensate the central value (asGG/n) = (376 MeV)^ and for 
the quark condensate (qq) eq. (|3.29|) . The 75-scheme stability is optimized at M = 
160 MeV. 
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of the residual scale dependence in the matching itself. We take the matching at 
0.8 GeV as the best compromise between the range of validity of the perturbative 
regime and that of chiral perturbation theory. 

We have also verified the consistency of the perturbative expansion in the short- 
distance regime by comparing leading order and the NLO results for Aq and A2 at 
fi = 0.8 GeV. The change in Aq turns out to be of about 10% and similarly for 
A2. These values confirm that the choice of the matching scale /i = 0.8 GeV is well 
within the perturbative regime of the QCD short-distance analysis. 

The values for Aq and A2 in table 4 should not be compared directly to the 
experimental values because they correspond to a specific choice for the quark and 
gluon condensates, as well as of M. A more complete discussion of the dependence 
of the results on the input parameters is presented in section 4.4. 



4.3 The Bi Factors 

For the purpose of comparison with the existing literature, we collect in table 5 
the Bi factors for the hadronic matrix elements of the operators in eq. ( p.2|) . The 
values of the Bi depend on the scale at which the matrix elements are evaluated, 
on the input parameters, on M, and on the 75-scheme; we have given in table 5 a 
representative example of their size. 

Only the first six operators are relevant for our present analysis and we postpone 
any comment on the electroweak penguin matrix elements when studying e'/e. 

The values of Bf^^ and B^^ show that the corresponding hadronic matrix ele- 
ments in the xQM are, once non-factorizable contributions and meson renormaliza- 

tion have been included, respectively twelve and three times larger than their VSA 

(2) (2) 

values. At the same time, B^ and B2 turn out to be about half of what found in 
the VSA. These features make it possible for the selection rule to be reproduced in 
the xQM. 



For comparison, in the 1/A^c approach of ref. [|10|, the inclusion of meson- loop 
renormalization through a cutoff regularization leads, at the scale of 1 GeV, to 

bJ°^ = 5.2, 5^°^ = 2.2 and = B^^'^ = 0.55 (see ref. H), a result that is not 

(2) (2) 

sufficient to reproduce the A/ = 1/2 rule. The value B\ = B2 = 0.55 in both 
the xQM (table 5) and the l/N^ approaches is remarkable, and yet a numerical 
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coincidence, since the suppression orginates from gluon condensate corrections in 
the xQM, whereas it is the effect of the meson loop renormahzation (regularized via 



explicit cut-off) in ref. |10 



The values of the penguin matrix elements (Q3) and {Q4) in the xQM lead to 
rather large Bi factors. In the case of Q3, the xQM result has the opposite sign of 
the VSA result and is negative. This is the effect of the large non-perturbative 
gluon correction. 

Regarding the gluon penguin operator Qe (and Q5), we find that the xQM 
gives a result consistent with the VSA (and the l/N^ approach), B^ {B5) being 
approximately equal to two for small values of the quark condensate and one at larger 
values. It is the quadratic dependence (to be contrasted to the linear dependence in 
the xQM) of the VSA matrix element for the penguin operators that it responsible 
for the different weight of these operators at different values of the quark condensate. 
The lattice estimate for these operators gives B^ = Bq = 1.0 ± 0.2 [^. 
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HV 


NDR 


u = 8 GeV 

IJj V/ • V.J V..' V 


ii = 1 GeV 

Laj J- • W V_J V 


it = 8 GeV 


u = 1 GeV 




10.6 


11.1 


10.6 


11.1 




2.8 


3.0 


2.8 


3.0 


B?^ 


0.52 


0.55 


0.52 


0.55 


B?^ 


0.52 


0.55 


0.52 


0.55 


B3 


-2.9 


-3.0 


-3.4 


-3.6 


B4 


1.8 


1.9 


1.3 


1.4 


B5 = Bq 


2.6^0.93 


2.7^0.98 


2.0^0.74 


2.2 0.78 


Bf^ 


3.2^2.3 


3.4 2.4 


3.0^2.2 


3.3^2.4 


^(0) 
-°8 


3.5^2.3 


3.8^2.5 


3.3^2.2 


3.6 2.4 


^(0) 
-°9 


3.9 


4.0 


3.6 


3.8 


-°10 


4.4 


4.7 


5.2 


5.5 


Bf^ 


2.6 1.4 


2.8^ 1.5 


2.6^ 1.4 


2.8^ 1.5 


^(2) 


2.1 ^ 1.4 


2.2^ 1.5 


2.0^ 1.4 


2.2 ^ 1.4 


^9 


0.52 


0.55 


0.52 


0.55 




0.52 


0.55 


0.52 


0.55 



Table 5: The Bi factors in the xQM (including meson-loop renormalizations) at two 
different scales: /i = 0.8 and 1.0 GeV. We have taken the gluon condensate at the 
central value of eq. ( |3.26|) , while the ranges given for B^_^ correspond to varying the 
quark condensate according to eq. ( p.27[ ). The results shown are given in the HV 
and NDR schemes for M = 180 MeV. 
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4.4 Dependence on the Input Parameters 



In this section we study the dependence of Aq and A2 on the quark and gluon 
condensates. Figs. 1, 2, 3 and 4 show such a dependence for Aq^j-, = 350 MeV and 
four representative values of M. 

Each figure contains the HV result (black lines) as well as the NDR one (grey 
lines). The spread between these two determinations is mostly due to A2 that 
contains an irreducible 75-scheme dependence, as explained in section 4.5. 

The gluon condensate dependence is represented by the three horizontal bands, 
the central one corresponding to the central value given in eq. ( p. 261) , while the other 



two bound the one standard deviation range. A larger gluon condensate leads to 
lower values of A2. 

The dependence on the quark condensate is represented by the points on a given 
line. The lenght of the lines corresponds to the range given in eq. (|3.27| ). 




Larger values of Aq are obtained for larger values of the quark condensate and 
smaller values of M. This is the effect of the gluon penguin operators Q^^e, whose 
matrix elements are proportional to the factor 

(4.5) 

in the HV scheme, and similarly (with 6 being replaced by 9) in the NDR scheme. 

The size of A2 is a sensitive function of the gluon condensate and the experimental 
value is approximately reproduced in both schemes by the central value of eq. ( p.26|) . 

Figs. 1, 2, 3 and 4 represent the main result of this paper. For any choice of M 
in the range 

M = 160 - 220 MeV (4.6) 

the experimental values of Aq and A2 are reproduced by our computation by different 
choice of the input parameters (qq) and (GG) that, however, remain within the 
ranges discussed in section 3.1. Values of M smaller than 160 MeV are still consistent 
with the experiments, while the suppression of the gluon penguin contributions as 
we increase M above 220 MeV would force us to quark condensate values outside 
the range of eq. ( p.27| ) in order to remain close to the experimental value for Aq. 



18 



3.5 
3 

> 25 

o 

^ 2 

oo 

< 

2 1.5 



X 
II 



1 

0.5 


-0.5 
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Figure 1: Dependence of Aq and A2 on {qq) and {GG) for Aq^j^ = 0.350 GeV, 
fx = 0.8 GeV, and M = 160 MeV. The black (grey) lines represent the HV (NDR) 
results by varying (qq) in the range of eq. ( |3.27| ) for fixed (asGG/n). The three 
bands correspond to varying (asGG/ir) in the range of eq. ( |3.26| ), with central lines 
correponding to the central value of (GG). The experimental values of Aq and A2 
are given by the cross hairs. The small dependence of A2 on the quark condensate 
is due to the the contribution of the electroweak penguins Qj g . 
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2 2.5 3 3.5 4 4.5 



A(I=0) X 10^7 (GeV) 

Figure 2: Same as in fig. 1 for M = 180 MeV. 



3.5 




2 2.5 3 3.5 4 4.5 



A(I=0) X 10^7 (GeV) 

Figure 3: Same as in fig. 1 for M = 200 MeV. 
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2 2.5 3 3.5 4 4.5 



A(I=0) X 10^7 (GeV) 

Figure 4: Same as in fig. 1 for M = 220 MeV. 
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4.5 75-SCHEME INDEPENDENCE 

The Wilson coefficients in eq. (|2.1|) depend at the NLO on the 75 scheme em- 
ployed |T^ , P^ . In the xQM also the hadronic matrix elements depend on the 
75-scheme used in the computation. The requirement that these two dependences 
balance each other to provide a 75-scheme independent result can be used in order 
to restrict the allowed values for the free parameter M. 

The scheme dependence is not strong in the amplitudes Aq and A2, and, there- 
fore, in this case, we can only restrict a range of values of M for which it is reasonably 
under control, but nothing more. The idea is however that when examining other 
observables with a stronger scheme dependence as, for instance, e' /e, we may find 
a substantial reduction of the scheme dependence for values of M that are still 
within the range determined by means of Aq and A2. A consistent picture for the 
whole of kaon physics would thus emerge, with different observables concurring in a 
consistent determination of the free parameter M as well as the range of the input 
parameters. 

In our estimate of the hadronic matrix elements, the 75-scheme independence of 
Aq turns out to be controlled by the gluonic penguins. In fact, the Wilson coefficients 
Z3_6 are systematically larger in the NDR scheme than in the HV scheme. On the 
other hand, in the NDR scheme the matrix elements of the gluonic penguins decrease 
with increasing M faster than in the HV scheme. As a consequence, it always exists 
a value of M for which the 75-scheme independence is achieved. 

Larger values of the quark condensate gives stability for larger values of M. From 
fig. 5 and 6 one can readily see that the difference between the two schemes remains 
below 20% for M between 140 and 240 MeV. 

The same is not true for the amplitude A2. In this case, the amplitude is con- 
trolled by the operators Qi and Q2 only so that there is no dependence on M and, 
accordingly, no intersection between the two schemes. The scheme dependence, 
however, stays well below 20%, as it can be seen in Fig. 7. 
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0.16 0.18 0.2 0.22 0.24 

M (GeV) 

Figure 5: The 75-scheme dependence of Ao is shown as a function of M. The black 
(grey) hne represent the HV (NDR) result. We use {asGG/ii) = (376 MeV)"^ and 
(qq) = —(220 MeV)^. For this value of the quark condensate the stability appear 
near M = 160 MeV. In the whole range of M shown, the 75-scheme dependence is 
always below 20%. 
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0.16 



0.18 0.2 
M (GeV) 



0.22 



0.24 



Figure 6: Same as in fig. 5 for {qq) = —(280 MeV)^. Tfie stability is moved at 
about M = 200 MeV. 




0.16 0.17 0.18 0.19 0.2 0.21 0.22 
M (GeV) 

Figure 7: Same as in fig. 5 for A2. 
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5 The Road to the A/ = 1/2 Selection Rule 

An instructive way of analyzing the relevance of the various contributions to the 
A/ = 1/2 rule is obtained by turning on in our computation each of them as we 
follow the historical steps that have lead to the present understanding of the rule 
(fig. 8). 

The point labeled by (1) in fig. 8 represents the theoretical prediction as obtained 
by considering the pure VSA matrix elements of Qi and Q2 without the short 
distance renormalization of the corresponding Wilson coefficients {fi = mw)- Point 
(2) represents the inclusion of the NLO renormalized Wilson coefficients, matched to 
the hadronic matrix elements at the scale fi = 0.8 GeV. This scale is large enough to 
make the renormalization-group analysis reliable while keeping the hadronic matrix 
elements in the chiral regime. 

As we can see, the value for Aq is far too small and that of A2 too large by a 
factor of two. 

The introduction of penguin operators (point (3)) goes in the direction of in- 
creasing Aq, but it leaves A2 unchanged. Their effect on Aq is not, at least in the 
xQM, as crucial as often claimed. 

The introduction of the electroweak penguins (Q7-10) little affects the CP- 
conserving amplitudes (point (4)), being suppressed by the smallness of their Wilson 
coefficients. Another isospin breaking contribution to the amplitude A2 comes from 
a long-distance effect, namely the mixing between 7r° and 77 particles. This contri- 
bution is evaluated to be 

Accordingly, we have a reduction of the amplitude A2 represented by point (5) which 
compensates for the effect of the electroweak penguins. Because of the smallness of 
Aq at point (5), the effect of ( |5.1|) is very small, though it reaches the 10% level in 
the final point (7). 

A crucial step toward the understanding of the selection rule is due to the (non- 
factorizable) gluon-condensate corrections (point (6)). They represent a genuine 
non-perturbative part of the computation. The isospin asymmetry generated by 
the electroweak operators is amplified in the right direction, improving dramatically 
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A2, that becomes close to its experimental value, while Aq is further increased with 
respect to point (5). The relevance of these contributions was first pointed out in 
ref. 0. 




Figure 8: The road to the AI = 1/2 rule: (1) Effect of the W- induced current- 
current matrix elements (Qi,2) with the neglect of short- distance QCD renormal- 
ization {fi = mw)] (2) {Qi,2) with the inclusion of the NLO Wilson coefficients 
at /i = 0.8 GeV. (3) inclusion of the gluon penguins (Qs-e); (4) inclusion of the 
electro-weak penguins {Qj-iq); (5) inclusion of the — rj mixing; (6) inclusion of 
gluon condensate corrections; (7) meson-loop renormalization. The results shown 
are those of the HV scheme with the central values (asGG/n) = (376 MeV)^ and 
{qq) = -(250 MeV)^ for a matching scale /i = 0.8 GeV and M = 180 MeV. The 
experimental values are given by the cross hairs. 

The meson-loop renormalization provides in our approach the final step toward 
the experimental results. The size of the relative renormalizations of Aq and A2 
goes in the right direction, being large for Aq and small for A2. The meson loops 
also introduce a renormalization scale dependence in the hadronic matrix elements 
to be matched with that of the Wilson coefficients. 
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As point (7) of fig. 8 shows, the AI =1/2 selection rule is well reproduced in 
the xQM, that provides not only values for the amplitudes Aq and A2 that are close 
to the experimental ones but also a satisfactory scale and 75-scheme independence 
of the estimate. Within 20% we also have scale independence in the matching range 
between 0.8 and 1 GeV. 

In drawing fig. 8 we have taken the gluon and quark condensates at the central 
values of eqs. ( p.26| )- p.271 ). Had we chosen, for instance, {asGG/n) = (372 MeV)"^ 
and (qq) = —(271 Me V)^ we would have exactly reproduced the experimental result. 
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Input Parameters 



parameter 


value 




0.9753 


Vus 


0.221 


sin^ Ow 


0.2247 


mz 


91.187 GeV 


mw 


80.22 GeV 


mi, 


4.8 GeV 


rric 


1.4 GeV 


fn fir+ 


92.4 MeV 


Ik = fK+ 


113 MeV 


rriTj- = (m7r+ + m^o)/2 


138 MeV 


rriK — ttiko 


498 MeV 


rrir, 


548 MeV 


Ax 




a(4) 


350 ± 100 MeV 


mu + md{i GeV) 


12 ± 2.5 MeV 


m 


-(200- 280 MeV)^ 


{(\,^GG / tt) 


(376 ± 17MeV)^ 



Table 6: Table of the numerical values of the input parameters. 
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